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EXECUTIVE  SUMMARY 


This  report  summarizes  work  we  have  undertaken  to  assemble  a  suite  of 
practical  mathematical  tools  that  allow  a  consistent  translation  of  standard 
weather  station  data  to  physiologically  relevant  input  streams  for  our  thermal 
strain  prediction  models.  Although  weather  station  measurements  of  air 
temperature,  relative  humidity/dew  point,  and  wind  speed  can  be  input  with 
minimal  pre-processing  into  most  thermal  strain  models,  the  solar  radiation 
component,  typically  measured  in  Watts/m^  by  an  up-looking  pyranometer,  can 
not.  This  report  draws  extensively  on  previously  published  work  and  describes 
processes  we  have  currently  implemented  to  use  global  pyranometer-based 
estimates  of  mean  radiant  temperature  as  the  common  solar  load  input  for  the 
Scenario  model,  the  USARIEM  heat  strain  model,  and  for  the  computation  of  the 
solar  radiation  sensitive  components  of  the  Wet  Bulb  Globe  Temperature 
(WBGT)  index. 

The  Mean  radiant  temperature  approach  averages  multiple  directional 
radiant  load  components  to  provide  a  single  uniform  surface  temperature  of  an 
imaginary  spherical  enclosure  surrounding  the  object.  For  estimating  mean 
radiant  temperature  from  weather  station  pyranometer  data  we  propose: 

Tmrk  =  [Tdbk^  +  (Py(1  +  agr)  +  S  gr  *  d  (Tg/  -  Tdbk"))/  2a  Ok 

where  Tmrk  is  the  mean  radiant  temperature  in  °K,  Py  is  the  Pyranometer 
reading  in  W*m'^,  Tdbk  is  dry  bulb  temperature  in  °K,  Tgrk  is  ground  surface 
temperature  in  °K,  8gr  is  ground  surface  emissivity  (0.7  to  0.9),  agr  is  ground 
surface  albedo  (  0.05  to  0.4  ),  and  a  is  the  Stefan  Boltzman  constant,  5.67  • 
10-® 

There  is  a  need  to  replace  the  category-based  solar  load  inputs  to  the 
USARIEM  heat  strain  model.  We  propose  a  continuous  function,  based  on 
mean  radiant  temperature  that  provides  a  smoothed  best  fit  to  the  categories 
across  the  natural  solar  radiation  domain: 


Rsoi  “  “0.071  •  (Tmr  “  Tdb  "*■  10.432  •  (Tmr  "  Tdb)  Watts 

where  Rsoi  is  the  total  radiant  load  impinging  on  the  outer  surface  of  clothing  in 
Watts,  Tmr  is  the  mean  radiant  temperature  in  °C,  and  Tdb  is  dry  bulb 
temperature  in  °C. 

Proposed  methods  for  computing  the  black  globe  and  natural  wet  bulb 
temperature  components  of  the  WBGT  index  from  mean  radiant  temperature  and 
the  other  standard  weather  measurements  are  described  and  documented  with 
appropriate  literature  citations. 
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INTRODUCTION 


OVERVIEW 

Strategies  for  quantifying  solar  load  are  incorporated  in  human  thermal  strain 
prediction  models  and  environmental  stress  indexes  such  as  the  wet  bulb  globe 
temperature  (WBGT)  index.  The  complex  asymmetrical  characteristics  of  most 
naturally  occurring  radiant  fields  has  led  to  the  application  of  a  range  of  sensor 
systems  as  potential  measurement  solutions  (2,  3, 4, 12, 13, 17).  Sensors  such 
as  the  black  globe  temperature  thermometer  and  more  elaborate  multi-directional 
radiation  sensing  systems  provide  the  basis  for  evaluating  solar  load.  Although 
a  simple  global  pyranometer  sensor  is  frequently  available  on  standard  military 
and  civil  weather  stations,  the  black  globe  thermometer  and  more  sophisticated 
directional  thermal  radiation  sensing  systems  are  not.  There  is  therefore  a 
continuing  need  to  assemble  practical  methods  to  translate  what  is  available 
from  standard  weather  stations  to  a  physiologically  relevant  context.  Simplified 
methods  to  produce  an  accepted  intermediate  metric  would  provide  a  common 
thread  that  facilitates  a  consistent,  coherent  computation  of  solar  load 
components  as  they  are  represented  both  in  the  indexes  and  the  more 
sophisticated  thermal  strain  prediction  models.  This  report  draws  extensively  on 
previously  published  work  and  describes  processes  we  have  currently 
implemented  to  use  global  pyranometer-based  estimates  of  mean  radiant 
temperature  as  the  common  solar  load  input  for  the  Scenario  model,  the 
USARIEM  heat  strain  model,  and  for  the  computation  of  the  solar  radiation 
sensitive  components  of  the  WBGT  index. 

Scenario  Model 

Mean  radiant  temperature  is  the  solar  load  input  parameter  for  the 
Scenario  model  (7,  8,  9).  Derivation  of  mean  radiant  temperature  from  the 
pyranometer  reading  completes  the  radiant  load  input  stream  requirements  for 
Scenario. 

USARIEM  Heat  Strain  Model 

Input  of  the  solar  load  component  of  heat  stress  for  the  original  USARIEM 
model  (15)  is  based  on  user  specification  of  one  of  four  broad  solar  categories: 
Darkness,  Cloudy,  Partly  Cloudy,  or  Full  Sun.  Internally,  the  category  is  then 
related  to  two  factors,  SIrF  and  CldF,  the  numerical  product  of  which  provides  a 
quantitatively  discreet  "band"  intended  to  represent  the  total  solar  radiation 
(Watts)  impinging  on  the  outer  clothing  surface.  This  value  is  subsequently 
multiplied  by  a  solar  efficiency  factor,  U,  that  accounts  for  the  effects  of  clothing 
and  wind  speed  on  solar  load  actually  received  by  the  body  (3).  The  category 
approach  is  a  reasonable  and  pragmatic  solution  to  situations  where  direct 
quantitative  measurements  of  the  radiant  components  are  unavailable  and 
reflects  the  philosophy  of  the  calculator-based  instantiation  of  the  model  which 
depended  entirely  on  externally  obtained  environmental  measurements  or 
approximations. 
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Recently,  the  USARIEM  model  has  been  implemented  in  a  hand  held  Heat 
Stress  Monitor  (HSM)  as  well  as  the  command  and  control  oriented  MERCURY 
test  bed  system  (10,11 ).  Both  of  these  systems  have  access  to  direct  or  derived 
measures^  of  ambient  solar  radiation.  In  order  to  exploit  these  data  and  automate 
the  model  input  process,  we  developed  schedules  that  allow  consistent 
selection  of  the  appropriate  solar  category.  Table  1  shows  the  categories  and 
the  schedule  by  which  they  were  initially  selected:  in  MERCURY  from  a 
pyranometer  measurement  of  global  solar  radiation,  and  in  the  HSM  from  mean 
radiant  temperature  (Tmr  )  derived  from  air  temperature  (Tdb),  wind  speed  (Vg) 
and  black  globe  (Tg)  measurements. 

Table  1.  Category-based  classification  scheduie  for  USARIEM  model  solar 
loads. 


Category 

SIrF 

(W) 

CldF 

SIrF’CldF 

(W) 

Pyranometer 

(W/m^) 

Tmr  ■  Tdb 
(°C) 

Darkness 

0.0 

4.0 

0 

<150 

<10 

3.5 

175 

150  to  399 

10  to  24.9 

Partly  Cloudy 

■SUMf 

3.0 

270 

400  to  700 

25  to  45 

Full  Sun 

150.0 

2.5 

375 

>700 

>45 

Statement  of  the  Problem 

It  was  found  during  sensitivity  studies  using  the  USARIEM  model 
implementation  for  the  Heat  Stress  Monitor  (HSM),  that  the  category  approach 
creates  disturbing  discontinuities  in  work/rest  cycle  and  maximum  safe  work  time 
outputs  when  the  solar  load  suddenly  shifts  to  a  new  category.  We  therefore 
decided  to  replace  the  category  based  solar  load  input  schedule  with  a 
continuous  function  that  would,  in  addition,  more  fully  exploit  the  currently 
available  solar  load-  related  measurements  available  to  both  the  HSM  and 
MERCURY'S  successor  system,  the  Operational  Medicine  Environmental  Grid 
Applications  (OMEGA)  Test  bed. 

OBJECTIVES  AND  APPROACH 

Objectives 

The  fundamental  objective  is  to  document  processes  we  have  used  to 
derive  mean  radiant  temperature  from  global  pyranometer  data  sets.  The  mean 
radiant  temperature  approach  averages  multiple  directional  radiant  load 
components  to  provide  a  single  uniform  surface  temperature  of  an  imaginary 
spherical  enclosure  surrounding  the  object.  A  key  objective  for  the  USARIEM 
heat  strain  model  applications  in  the  HSM  and  OMEGA  is  to  document 
methodology  for  an  empirically  derived  continuous  solar  load  input  parameter, 
Rsoi,  that  is  based  on  mean  radiant  temperature,  Tmr,  and  essentially  overlays 
that  model's  category-based  solar  load  inputs.  The  intent  is  to  preserve,  for 
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better  or  worse,  the  quantitative  trends  explicit  in  the  category-based  handling 
of  solar  load  and  to  achieve  a  historically  consistent  but  significantly  smoother 
quantitative  representation  of  this  model  input  across  the  solar  load  domain.  The 
primary  objective  for  the  WBGT  index  components  is  to  document  those 
mathematical  relationships  we  have  adopted  for  estimating  black  globe  and 
natural  wet  bulb  temperatures  using  Tmr  • 

Approach 

The  approach  is  to  employ  regression  analysis  to  establish  a  smooth  fit  to 
the  category  products  using  Tmr  as  the  common  intermediate  parameter  for 
solar  load  inputs  to  both  the  HSM  and  OMEGA  model  instantiations.  Because 
Tmr  is  easily  derived  from  the  HSM's  black  globe,  air  temperature  and  wind  speed 
measurements  using  standard  published  equations  (1),  this  report  will  focus  on 
methodologies  for  approximating  Tmr  from  pyranometer  data.  A  significant 
collateral  benefit  of  pyranometer-derived  Tmr  is  that,  in  conjunction  with  the  other 
standard  meteorological  measurements  of  air  temperature,  humidity,  and  wind 
speed,  it  is  possible  to  estimate  black  globe  temperature  (Tg  ),  and  the  naturally 
convected  wet  bulb  temperature  (Tnwb  )  which  are  the  two  primary  components  of 
the  WBGT  index.  In  spite  of  known  limitations,  WBGT  continues  to  be  the 
preferred  heat  stress  management  tool  within  the  military  preventive  medicine 
community.  Its  inclusion  as  an  ancillary  product  in  more  rational  predictive  model 
output  streams  (14)  is  desirable  because  the  ability  to  demonstrate  inadequacies 
of  this  culturally  entrenched  frame  of  reference  is  essential  to  wider  acceptance 
of  rational  models.  Our  approach,  based  on  pyranometer-derived  Tmr ,  is  thus 
intended  to  provide  a  coherent  sensor  fusion  methodology  for  linkage  of 
standard  meteorological  sensor  systems  to  the  WBGT  and  to  the  more  rigorous 
thermal  strain  prediction  models. 


METHODS 

ANALYTICAL  PROCEDURES  AND  SENSOR  SYSTEMS 

This  work  is  intended  to  document  the  historical  and  analytical  basis  for 
changes  in  the  solar  load  related  software  modules  currently  or  soon  to  be 
instantiated  in  USARIEM's  thermal  strain  predictive  models.  The  methodologies 
developed  in  the  results  section  are  based  on  analysis  of  data  obtainable  from 
sensor  systems  commonly  available  on  weather  station  platforms:  Air 
temperature,  relative  humidity,  wind  speed,  and  solar  radiation  sensors 
comprise  the  core  suite  of  sensors  for  current  thermal  strain  prediction  models. 
The  analytical  focus  is  the  extraction  of  physiologically-relevant  solar  load 
information  from  pyranometer  data.  In  clear  sky  situations,  readings  from  these 
sensors  range  typically  from  slightly  below  0  at  night  to  more  than  1000  W*m'^ 
at  mid  day.  Data  from  up-looking  pyranometers  are  available  from  automated 
weather  station  systems  such  as  the  Remote  Automated  Weather  Stations 
(RAWS)  system  at  Eglin  AFB,  Florida  and  USARIEM's  new  portable  RAWS 
system  that  supports  OMEGA  and  the  Warfighter  Physiological  Status 
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Monitoring  (WPSM).  Statistical  evaluations  and  regression  analyses  required  for 
this  effort  were  accomplished  using  tools  available  in  Microsoft  Excel/Office  2000 

TM 


RESULTS 

USARIEM  MODEL  AND  PYRANOMETER  DATA 

A  direct  relationship  between  the  pyranometer  reading  and  the  SlrF«CldF 
product  was  derived  using  a  polynomial  regression  fit  as  shown  in  Figure  1 .  As  a 
smoothed  representation  of  the  SlrF*CldF  product,  the  variable  Rsoi  quantifies  the 
extra  radiation  impinging  on  the  outer  surface  of  the  clothing  in  Watts,  for 
conditions  where  air  and  radiant  temperature  differ.  It  is  important  to  note  that 
Rsoi  reflects  significant  assumptions  about  the  effects  of  high  sun  angles  on 
projected  area  of  a  standing  person. 

Figure  1 .  Rsoi :  A  polynomial  fit  of  category-based  solar  load  to 
pyranometer  readings. 


MEAN  RADIANT  TEMPERATURE 

The  pyranometer  and  air  temperature  measurements  from  the  Remote 
Automated  Weather  Stations  (RAWS)  are  used  as  the  basis  for  estimating  Tmr- 
Using  the  pyranometer  reading  to  compute  a  mean  radiant  temperature  provides 
significant  advantages  over  the  direct  empirical  fit  described  above.  It  allows 
consideration  of  other  indirect  but  related  radiant  load  components  and  has  more 
general  utility  as  an  input  parameter  for  predictive  models  such  as  Scenario  and 
for  the  WBGT  derivation  algorithms.  Because  it  is  unrealistic  to  assume  that  all 
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of  the  radiant  load  is  coming  from  the  sky  hemisphere  and  is  measured  by  the 
up-looking  pyranometer,  additional  measurements  or  reasonable  assumptions 
are  in  order  to  account  for  reflected  solar  radiation  and  upwelling  long  wave 
radiation  from  the  heated  ground  surface.  An  isotropic  mean  radiant  field  is 
estimated,  as  the  sum  of  the  4%  contributions  of  each  of  these  elements.  Note 
that  temperature  units  for  dry  bulb  temperature,  ground  temperature,  and  mean 
radiant  temperature  calculations  are  degrees  Kelvin  (°K  =  +  273.15). 

Figure  2  is  a  schematic  of  a  mean  radiant  temperature  (Tmrk) 
representation  of  the  radiation  received  by  an  above  ground  radiantly  black 
object  at  air  temperature  (Tdbk)-  The  object  receives  direct  and  ground  reflected 
solar  radiation  and  long  wave  radiation  from  the  ground. 

Figure  2.  Schematic  of  radiation  received  by  an  above  ground  object. 


The  angle  of  incidence  of  the  direct  solar  beam  radiation  is  0.  The  direct  (Id)  and 
diffuse  (Idif)  solar  radiation  are  incident  on  the  projected  up-looking  horizontal 
area  (Aup)  of  the  object.  The  total  intensity  received  is  (Id’cos  0  +  Idif)  which  is 
assumed  to  equal  the  incident  radiation  intensity  (Py)  measured  by  a  skyward 
looking  horizontal  pyrometer. 

The  solar  radiation  reflected  from  the  ground  to  the  underside  of  the  object  is 
Adown'Ir,  where  Ir  is  the  direct  and  diffuse  solar  radiation  reflected  by  the  ground. 
It  is  represented  as  agr*Py  where  agr  is  the  ground’s  albedo  or  its  hemispherical 
solar  reflectance.  Representative  values  for  the  albedo  of  some  common 
surfaces,  from  Holman  (6),  are  given  in  Table  2. 
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Table  2.  Albedos  for  some  natural  surfaces  (From  Holman,  1997). 


Surface 

Albedo 

Water 

0.3-0. 4 

Black  dry  soil 

0.14 

Black  moist  soil 

0.08 

Gray  dry  soil 

0.25-0.30 

Gray  moist  soil 

0.10-0.12 

Desert  loam 

0.29-0.31 

Bright  fine  sand 

0.37 

Snow 

0.4-0.85 

Green  grass 

0.26 

Further,  the  underside  of  the  reference  object  may  also  receive  radiant  heat  (qgr) 
from  the  ground. 


qgr  =  8gr-a(Tgrk"-Tdbk")  Watts  (Eq.1) 

where  Sg,  is  the  long  wave  (low  temperature)  emissivity  of  the  ground.  For 
opaque  materials  absorptivity  equals  emissivity  and  some  representative  values 
from  Holman  (6)  are  listed  in  Table  3. 

Table  3.  Absorptivities  /  emissivities  of  various  surfaces  to  solar  and  long 
wave  radiation  ( From  Holman,  1997). 


Surface 

solar 

long  wave 

Aluminum,  highly  polished 

0.15 

0.04 

Copper,  highly  polished 

0.18 

0.03 

Tarnished 

0.65 

0.75 

Cast  iron 

0.95 

0.21 

Stainless  steel,  no.  301,  polished 

0.37 

0.60 

White  marble 

0.46 

0.95 

Asphalt 

0.90 

Red  brick 

0.75 

0.93 

Gravel 

0.29 

0.85 

White  paint 

The  radiant  energy  depicted  in  Figure  2  can  be  mathematically  represented  as; 

(Aup+Adown)  *  O  (Tmrk  ~Tcibk  )  ~  Ayp  (Id  ‘COS  0  +  Idif)  Adown  (Ir  "'"Qgr)  (EQ'2) 

The  projected  areas  facing  up  and  down  must  be  equal  (Aup  =  Adown  =  A) 
simplifying  Eq.  2  to: 
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2A*<J  (Tmrk^  —  Tdbk^)  “  A(ld  ’COS  6  +  Icjif)  +  A  (If  ^Clgr) 

Then  substituting  the  Pyrometer  reading  (Py)  and  simplifying: 

'  2cr  (Tmrk  ~  T(jbk  )  ”  Py  (1  "*■  3gr)  +  £  gr  *<7  (Tgrk  ~  Tdbk  ) 

or 

Tmrk  “  [Tdbk  (Py(1'*’  ^gr)  £gr*cr  (Tgrk  ~  Tdbk  ))/  2cT  ] 

where: 

Tmrk  =  mean  radiant  temperature,  °K 
Py  =  pyranometer  reading,  W»m'^ 

Tdbk  =  dry  bulb  temperature,  °K 
Tgrk  =  ground  surface  temperature,  °K 
£gr  =  ground  surface  emissivity 
a-gr  =  ground  surface  albedo 
a  =  Stefan  Boltzman  constant,  5.67  x  1 0'®  W*m'^*  °K'^ 

Figure  3.  Effects  of  ground  temperature  on  estimated  mean  radiant 
temperature  at  constant  30  °  air  temperature. 
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°K  (Eq.5) 


Ground  surface  temperature  measurements  are  provided  by  RAWS  systems  and 
reasonable  approximations  for  bare  ground  albedo  and  long  wave  emissivity  can 


be  obtained  from  Tables  2  and  3.  Note  that  if  the  ground  albedo  is  assumed  to 
be  0  and  the  ground  temperature  equals  air  temperature  eq.  5  simplifies  to: 

Tmrk  =  (Tdbk^  +  Py)/  2a  Ok  (Eq.6) 


USARIEM  MODEL  AND  MEAN  RADIANT  TEMPERATURE 

Using  Tmr  -  Tdb  categories  from  Table  1,  regression  analysis  yields  a 
continuous,  smoothed  function  that  overlays  the  USARIEM  model's  category- 
based  solar  load  as  shown  in  Figure  4 . 

Figure  4.  Rsoi:  A  polynomial  fit  of  category-based  solar  load  to  Tmr  -  Tdb  • 


Tmr-Tdb  (°C) 


Polynomial  regression  yields: 

Rsoi  =  -0.071  (T,nr  -  Tdb  f  +  10.432  (T^r  -  Tdb)  Watts  (Eq.7) 

The  continuous  variable  Rsoi  here  represents  a  "plug-in"  replacement  for  the 
SIrF’CldF  product  in  the  USARIEM  heat  strain  model  algorithm  that,  along  with 
the  solar  efficiency  factor,  U,  is  used  to  compute  the  actual  radiant  load  on  the 
body.  While  it  appears  to  be  a  reasonable  approximation  of  the  USARIEM 
model's  handling  of  solar  load,  it  is  probably  appropriate  to  characterize  the 
magnitude  of  departures  from  a  more  generalized  rational  estimate  such  as 
Gagge's  Effective  Radiant  Flux  (ERF). 
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Effective  Radiant  Flux  (ERF) 

Effective  Radiant  Flux  (Santee  and  Gonzalez,  1988)  is  similar  to  Rsoi  in 
that  it  represents  solar  load  as  a  function  of  air  temperature  and  mean  radiant 
temperature,  but  it  also  includes  terms  for  clothing  absorptivity  and  radiating 
surface  ar;ea.  We  compute  whole  body  ERF  as  follows: 

ERF=  a‘a*Ar/AD  •(Tmrk'^-Tdbk^-AD  Watts  (Eq.8) 

Where, 

Tdbk  =  dry  bulb  temperature,  °K 

a  =  absorptivity  of  clothing,  here  0.9 

cr  =  Stefan  Boltzman  constant,  5.67  x  10'®  W*m'^*°K'^ 

Tmrk  =  mean  radiant  temperature,  °K 
Ad  =  Dubois  body  surface  area,  here  1 .8  m^ 

Ar/  Ad  =  ratio  of  radiating  surface  to  total  body  surface  area,  0.72 

Figure  5  shows  the  relationship  between  ERF  and  the  SlrF*CldF  product. 

Figure  5.  Comparison  of  Effective  Radiant  Flux  (ERF)  and  Rsoi  as  a 
function  of  Tmr  -  Tdb  • 


Aside  from  the  fact  that  clothing  effects  on  net  solar  load  burden  in  the  USARIEM 
model  are  computed  following  the  estimate  of  Rsoi,  it  is  clear  that  significant 
implicit  assumptions  about  view  angle,  solar  altitude,  and  posture  are  aggregated 
in  the  category  factors  The  difference  between  the  computed  ERF  and  Rsoi  Is 
shown  in  Figure  6. 
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Figure  7  depicts  the  polynomial  regression  function  that  yields  an  adjustment 
factor/multiplier  for  ERF  that  quantifies  the  difference  between  ERF  and  Rsol 
across  the  solar  load  domain.  The  ERF  multiplier  ranges  from  1 .4  to  roughly 
0.45  at  the  highest  solar  load. 

The  expression  modifying  ERF  and  relating  ERF  to  Rsoi  thus  becomes: 

Modified  ERF  = 

Rsoi=  ERF  •  (0.00005  (T,T,r-Tdb)^- 0.0161  (Tmr-Tdb)  +  1.4)  Watts  (Eq.9) 
Figure  8  shows  the  modified  ERF  across  the  solar  load  domain. 


Figure  8.  Modified  ERF  as  a  function  of  Tmr  -  Tdb  • 


Although  the  preceding  exercise  characterizes  the  relationship  between  ERF  and 
Rsol,  and  would  allow  direct  compatibility  with  the  USARIEM  heat  strain  model 
input  stream,  applications  in  OMEGA  and  the  Heat  Stress  Monitor  will  use  the 
exclusively  empirical  fit  described  previously,  Eq.  7: 

Rsol  =  -0.071  (Tmr  -  Tdb  ?  +  10.432  (Tmr  -  Tdb)  Watts 


COMPUTATION  OF  THE  WBGT  COMPONENTS 

The  addition  of  a  Wet  Bulb  Globe  Temperature  (WBGT)  capability  for 
OMEGA  does  not  represent  a  cutting  edge  technological  improvement  for 
thermal  injury  risk  assessment.  It  does,  however,  provide  a  useful  resource  to 
support  comparative  and  analytical  evaluations  of  this  widely  used  heat  stress 
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index.  It  also  provides  an  important  frame  of  reference  and  credible  historical 
thread  for  potential  OMEGA  users  who  are  currently  committed  to  WBGT-based 
heat  stress  management  doctrine  (16). 


Part  of  the  appeal  and  resilience  of  WBGT  methodology  is  the  simplicity  of  the 
basic  measurements  and  the  computation  of  the  index.  The  WBGT  index  for 
outdoor  environments  is  computed  as  follows: 

WBGT=0.7Tnwb  +0.2Tg  +0.1Tdb  °Cor°F  (Eq.10) 


Where; 

Tnwb  =  Natural  Wet  Bulb  Temperature  °CorT 
Tg  =6  inch  Black  Globe  Temperature  °C  or  °F 
Tdb  =  Air  Temperature  °C  or  °F 


Because  two  components  of  the  WBGT,  naturally  convected  wet  bulb 
temperature,  Tnwb  >  and  black  globe  temperature,  Tg ,  are  not  measured  by 
standard  weather  stations,  these  values  must  be  derived  from  the  available 
measured  data.  The  original  conceptual  advantage  of  WBGT  in  terms  of 
simplicity  is  certainly  lost  when  these  components  must  be  derived  from  standard 
measures  using  iterative  computational  methods,  but  the  processes  can  be 
automated  in  software.  It  should  be  noted  that  recent  development  of  a 
regression  based  equation  that  uses  weather  station  data  to  compute  an 
Environmental  Strain  Index  (ESI)  has  demonstrated  very  high  correlation  with 
measured  WBGT  values  (12).  This  may  ultimately  prove  to  be  an  acceptable 
substitute  for  the  more  computationally  demanding  solutions  described  below. 


Natural  Wet  Bulb  Temperature 

In  the  current  implementation  of  OMEGA,  the  naturally  convected  wet  bulb 
temperature  is  computed  using  an  iterative  process  to  solve  a  previously 
described  equilibrium  state  equation  (5).  A  first  guess  for  Tnwbk  (Tnwbk-iEsi , 
starting  typically  at  Tdbk)  is  substituted  into  equation  1 1  and  both  sides  evaluated. 
Tnwbk-TEST  IS  incrementally  decremented  until  both  sides  converge  on  a  value  that 
satisfies  the  equality,  and  Tnwbk-iEST  equals  the  true  Tnwbk  •  The  process  uses  the 
derived  Mean  Radiant  Temperature  (Tmrk),  and  the  air  temperature  (Tdbk),  wind 
speed  (Va ),  and  relative  humidity  (RH)  measurements.  From  Gonzalez  et  al.: 

hc(Tdbk"Tnw  bk)  Bnwb  *  (Tmrk  "Tpwbk  )“  he(Psnwb  "  (RH/100)*Psdb)  (Eq.11) 


Where: 

he  =  convective  heat  transfer  coefficient,  W*m’^*  °K'^ 

Tdbk  =  dry  bulb  temperature,  °K 

Tnwbk  =  natural  wet  bulb  temperature,  °K 

Snwb  =  emissivity  of  the  surface  of  the  wetted  wick,  ~  0.5 

a  =  Stefan  Boltzman  constant,  5.67  x  10  ®  W'm‘^*°K'^ 

Tmrk  =  mean  radiant  temperature,  °K 
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he  =  evaporative  heat  transfer  coefficient,  WTn"^«Torr^ 

Psnwb  =  saturation  water  vapor  pressure  atTnwb,  Torr 
RH  =  %RH 

Psdb  =  saturation  water  vapor  pressure  at  Tdb,  Torr 

The  convective  heat  transfer  coefficient  is  computed  from  wind  velocity,  Va  ,  in 
m*s-'': 


he  =  42.024  W*m-2‘°K-''  (Eq.12) 

And  the  evaporative  heat  transfer  coefficient  is  related  to  the  convective  heat 
transfer  coefficient  by  the  Lewis  relation: 


he  =2.2  he  W*m'^-Torr^  (Eq.13) 

Saturation  water  vapor  pressures  at  Tnwb  and  Tdb  are  computed  using  a  base  10 
version  of  Antoine's  equation  for  temperatures  in  °C  (Santee  and  Gonzalez, 
1988). 

.|q(8.10765- 1750.286/(7+235))  (EO.M) 


Black  Globe  Temperature 

In  the  current  implementation  of  OMEGA,  black  globe  temperature  is 
computed  using  a  previously  described  equilibrium  state  equation  for  the  globe, 
where  radiative  heat  gain  equals  convective  heat  loss  (ASHRAE,  1984): 


or 


hrg  (Tmrk  “  Tgk)  —  heg  (Tgk  —  Tdbk) 

Tgk  =  (heg  ‘Tdbk  hrg  ’Tirirk)  /  (heg  hrg) 


(Eq.15) 


Where: 

Tgk  =  black  globe  temperature,  °K 
hrg  =  linear  radiation  transfer  coefficient  for  globe 
=  Sg  •  (J  (Tmrk  +  Tdbk  )  •  (Tmrk^  +  Tdbk^  ),  W»m‘^  • 
heg  =  convective  heat  transfer  coefficient 
=  6.32  D-°-^  •  Va°-®  ,  W‘m-2  • 


8g  =  emissivity  of  black  globe,  -0.95 
a  =  Stefan  Boltzman  constant,  5.67  x  10‘®  W+m'^+^K"^ 
Tdbk  =  dry  bulb  temperature,  °K 
Tmrk  =  mean  radiant  temperature,  °K 
D  =  diameter  of  globe,  m 
Va  =  air  velocity,  m-s'^ 
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DISCUSSION 


MODEL  INTEGRATION 

As  the  federation  of  warfighter  predictive/performance  models  required  to 
support  Army  science  and  technology  objectives  becomes  more  tightly 
integrated,  there  is  a  growing  need  to  define  and  document  practical  methods  for 
handling  environmental  measurements.  The  ability  to  leverage  existing 
meteorological  resources  as  input  streams  for  model  test  and  evaluation  expands 
research  opportunities,  allows  early  identification  of  meteorological  data  gaps  or 
sensor  inadequacies,  and  ultimately  contributes  to  the  effective  definition  of 
meteorological  requirements  for  highly  integrated  operational  systems  such  as 
Land  Warrior. 


EXPLOITING  EXISTING  INFRASTRUCTURE 

Our  efforts  here  focus  on  the  solar  radiation  component  of  the  natural 
environment  and  the  derivation  of  a  commonly  accepted,  physiologically  relevant 
metric  that  would  support  a  range  of  indexes  and  thermal  strain  prediction 
models.  The  methodology  proposed  for  derivation  of  Tmr  from  pyranometer  data 
represents  a  considerable  compromise  from  the  most  robust  biophysical 
evaluations.  Nevertheless,  it  does  provide  a  practical  solution  to  a  data  gap  that 
is  not  likely  to  be  otherwise  resolved  in  the  near  future.  Pyranometers,  based  on 
silicon  photodiode  technology,  are  widely  used  solar  radiation  sensors.  Data 
from  up-looking  pyranometers  are  available  from  automated  weather  station 
systems  such  as  the  Remote  Automated  Weather  Stations  (RAWS)  system  at 
Eglin  AFB,  Florida  and  USARIEM’s  new  portable  RAWS  system  that  supports 
OMEGA  and  the  Warfighter  Physiological  Status  Monitoring  (WPSM)  project. 


CONCLUSIONS 


There  is  a  growing  need  to  exploit  existing  standard  weather  station  data 
resources  to  support  emerging  requirements  for  real  time-modeling  of  Warfighter 
thermal,  hydration,  and  performance  status. 

Combined  with  standard  measurements  of  air  temperature,  humidity,  and  wind 
speed,  a  practical  method  to  derive  mean  radiant  temperature  from  weather 
station  pyranometer  readings  provides  a  basic  input  stream  for  thermal  strain 
prediction  models  and  provides  opportunities  for  physiologically  relevant  sensor 
fusion. 

The  proposed  expression  relating  pyranometer  readings  to  mean  radiant 
temperature  provides  a  common  thread  for  consistent  utilization  of  standard 
meteorological  data  for  automated  input  of  solar  load  to  the  Scenario  and 
USARIEM  heat  strain  models. 
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The  proposed  expression  relating  pyranometer  readings  to  mean  radiant 
temperature  provides  a  common  thread  for  consistent  utilization  of  standard 
meteorological  data  for  rational  estimation  of  the  individual  components  of  the 
traditional  Wet  Bulb  Globe  Temperature  (WBGT)  index. 

The  proposed  continuous  function  relating  mean  radiant  temperature  to  the 
category-based  solar  load  input  in  the  USARIEM  heat  strain  model  provides  a 
practical  solution  for  linkage  of  that  model  to  standard  weather  data  resources. 
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